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SUMMARY 

Stresses and deformations in thick, short, composite cylinder subjected to 
axially varying radial pressure are studied. Effect of slippage at the inter- 
face is examined. In the NASTRAN finite element model, multipoint constraint 
feature is utilized. Results are compared with theoretical analysis and SAP-IV 
computer code. Results from NASTRAN computer code are in good agreement with 
the analytical solutions. Results suggest a considerable influence of inter- 
facial slippage on the axial bending stresses in the cylinder. 


INTRODUCTION 


The analysis of composite (multilayered) cylinders and other axisymmetric 
bodies with shrink (or press) fit forms an extremely useful class of problem 
in the design of pressure vessels and various machine components. For example. 
Figure 1 illustrates a special slipping pipe anchor (Ref. 1) consisting of the 
clamped (or press fit) support. Figure 2 illustrates a typical bar drawing die 
composed of an outer hardened steel cylinder press fitted on the inner carbide 
insert. The loading, in both cases, on the inside surface of the cylinder is 
variable axial pressure dependent on the seismic motions (Fig. 1) and percent- 
age reductions (Fig. 2) respectively. This paper attempts to present the prob- 
lem in a generalized way and indicates, in quantitative terms, the effect of 
interfacial slippage on the various stress components. 


SYMBOLS 


E Young's modulus of elasticity, MPa 

V Poisson's ratio (Note: Suffix 1 refers to inner cylinder and 2 refers 

to outer cylinder) 

u^ Radial displacements, cm. 

u^ Axial displacements, cm. 

L Length of the cylinder, cm. 
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Z Axial coordinate 

y,r Radial coordinate 

P Internal pressure, MPa 

k Scale factor for pressure 



Axial bending stress, MPa 
Hoop stress, MPa 
Radial stress , MPa 
Inside radius, cm. 
Interface radius, cm. 
Outside radius, cm. 
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ANALYTICAL MODELS 


Figure 3 shows the NASTRAN (Ref, 2) finite element model of the composite 
cylinder. The inner cylinder is of carbide and the outer cylinder is of steel. 
The model consists of 210 nodes and 168 axisymmetric trapezoidal ring (CTRAPRG) 
elements. Dimensions and material properties are typical of a bar drawing die 
used in actual steel drawing process. 

The loading consists of an axially varying radial pressure defined by the 
equation P « k(Lz - z2). This pressure was input as discrete load on the in- 
side surface using the ’FORCE* card of NASTRAN, 

Nodes 91 through 120 are nodes along the interface of the carbide and 
steel. This interface needs special consideration in the finite element model. 
In practice, the case is shrinkfitted on the carbide insert, and a large amount 
of compressive preloading is generated in the relatively brittle carbide insert. 
This compressive preloading reduces the probability of cracking of the insert 
during operation. 

Since there is no welding (or permanent connection) along the interface, 
during nonuniform loading, slippage will occur along the mating surfaces. At 
the interface, compatibility of the radial displacement provides: 


u \ = u \ for 0 ^ Z £ L 

/ carbide / steel 

The axial displacement could vary depending on the amount of slippage. 

Thus 


u 

z 


) 


carbide 




) Steel 


NASTRAN (Ref, 3) has a very useful option available through the MPC command. 
It defines the multipoint constraint of the form: 


2A u 
j J ^ 


0 where Aj’s are the real coefficients and u. 
placements. 


are the nodal dis- 


using this option, the present analysis was carried out on NASTRAN level 
15,5 through the UCC 1108 System. 


Figure 4 shows the finite element model for the SAP-IV computer code 
(Ref. 4), This model is similar to the NASTRAN model - however, no slippage 
can be allowed along the interface. Due to this, the compatibility of defor- 
mation along the interface does not represent the actual deformation mode. 


THEORETICAL ANALYSIS 

A theoretical analysis of this problem is carried out by the author*. In 
this analysis, the governing equations of three dimension elasticity are 


*Being published as a separate paper 
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satisfied on the curved surfaces using the series approach of Lee (Ref S') 
The residual elasticity problem is then solved to partially satisfy th; * 

ends of the cylinder. The compatibility of dis- 
it! all ^ H a t e interface is included in such a manner that axial slip 
as follow^: deformation. The resultant stress components derived are 
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RESULTS AND DISCUSSIONS 


Fig. 5 shows the variation of axial stresses ( 022 ^ along the 

length of the inner cylinder. The axial stresses reach their maximum value at 
the center of the cylinder and the variation is symmetrical about Z = L/2 due 
to the symmetry of the loading. The results from tl>e NASTRAN model are in 
good agreement with the analytical solution. Further, it may be noted that 
t e results are in better agreement on the inside surface than the outside sur- 
face of inner cylinder - the maximum deviations are 2.8% and 8% respectively. 
This is due to the fact that the inner surface is farther away from the slip- 
ping interface thereby less sensitive to the boundary effects. 
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Variations of the hoop stresses (Figure 6), radial stresses (Fig. 7) and 
deflections (Fig. 8-9) also indicate good agreement between the analysis and 
NASTRAN model. At the Interface, in Fig. 7, the radial stresses on the out- 
side surface of carbide and inner surface of steel are identical from the 
analytical solution. In Figure 9, the length is twice the value given in 
Fig. 3 while the remaining dimensions are the same as Fig. 3. 

The effect of slippage at the interface, in quantitative terms, is demon- 
strated by the values of stresses from SAP IV, NASTRAN and the analytical solu- 
tion (Table 1). As is indicated by this table, the effect is more pronounced in 
the values of a - since the slippage makes the bending process more symmet- 
rical. 


CONCLUDING REMARKS 


The stresses and deformations in composite cylinder subjected to axially 
varying radial pressure are examined using the NASTRAN finite element models. 

The MPC feature of NASTRAN is very useful in simulating accurately the slipping 
interface between carbide and steel. Results from the NASTRAN finite element 
models and the analytical solutions are found to be in good agreement. 

TABLE 1 


Stresses in Carbide (Inner) Cylinder (Units: MPa) 


Description 

Analytical 

Solution 

NASTRAN 
F.E. Model 

SAP IV 
F.E. Model 

At radius = 5. 08 cm. 

a 

-124.6 

-121*0 

-136*2 

j zz 

Q 

CD 

CD 

362.3 

366*2 

355*0 

a 

-223*4 

-221*8 

-220*6 

rr 

At radius = 6*99 cm. 

a 

zz 

110*93 

102*0 

26*2 (Regular) 
87*5 (Modified*) 

®90 

279.7 

275*2 

279*2 

a 

-69*9 

-74*4 

-72*4 

1 rr 

(^Modified using, in steel 

at interface. 

one layer of a 

material having low 

modulus of elasticity in z 

direction) 


165 


REFERENCES 


1 . 

2 . 

3. 

4. 


Patel, Y A., Cho, F. L.; Diinopoulos, A. P.: Slipping Pipe Anchor - A Novel 

Piping in High Seismic Intensity Zones - ASME Paper 75-PVP-58 
ASHE Conf. on Pressure Vessels and Piping, June 1975. 

C.W. McCormick (Edit.) - The NASTRAN User’s Manual - NASA SP-222(01) - 
National Aeronautics and Space Administration, June 1972. 

Lee, C. W.: A Theory of Thick Walled Cylinder Under Axisymmetrlc Loading - 

roceedings of the Fourth U. S. National Congress of Applied Mechanics 
(ASME) N.Y., 1962, p. 662. «PP-txea necnanics, 

Bath, K. J.; Wilson, E. L. and Peterson, F. E.: SAP-IV - A Structural 

Program for Static and Dynamic Response of Linear Systems. 

EERC Report 73-11, University of Calif., Berkeley, June 1973. 


166 


Insulation 



167 


Fig. 1 Slipping pipe anchor. 
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Section of a typical draw die and workpiece. 
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Fig. 3 Nastran finite element model of the composite c^inder. 
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Fig. 6 Variation of hoop stresses with axial distance. 




Fig. 7 Variation of radial stresses with axial distance. 
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— ^-nNastran Model 
Analytical 





Fig. 9 Variation of radial deflection. 
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